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ABSTRACT 


The  slow  pyrolysis  of  lithium  aluminum  hydride  proceeds  by  a 
successive  dehydrogenation  of  the  fuel  in  three  stages  over  a  temperature 
range  between  210°  and  460°C.  The  combustion  of  lithium  aluminum 
hydride  follows  a  simular  multistage  process  in  which  the  pyrolytic 
dehydrogenation  is  followed  by  an  afterburning  of  the  metal  substrate. 
Experimental  studies  were  performed  to  determine  the  delay  times 
between  successive  stages  for  solid  mixtures  of  lithium  aluminum 
hydride  with  selected  oxidizers.  Time  resolved  spectrographic  tech¬ 
niques  were  used  to  determine  the  duration  of  these  stages.  Studies 
were  performed  using  Li202  and  NH4N03  oxidizers  over  a  particle  size 
range  between  75-300  microns  and  O/F  mixture  ratios  of  between  0.  5 
and  3.  0.  Ignition  delay  times  for  the  initial  reaction  were  found  to 
vary  between  5  and  40  milliseconds  for  the  first  stage  and  10  to 
120  milliseconds  for  the  second  stages.  An  analytical  model  was  devel¬ 
oped  to  describe  the  observed  multistage  combustion  process.  The 
model  and  resulting  PARTIG  computer  program  utilize  the  nodal 
approach  with  stepwise  solution  of  the  equations  for  heat  and  mass 
transfer.  The  analytically  predicted  ignition  delay  times  for  a  50p 
particle  agreed  favorably  with  the  experimentally  obtained  times. 
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I.  INTRODUCTION 

Advanced  rocket  propellant  combinations  include  the  use  of  light 
metal  hydride  fuels.  The  physical  and  chemical  mechanisms  important 
to  the  combustion  of  metal  hydrides  differ  substantially  from  those 
associated  with  the  combustion  of  nacent  metals.  For  metal  hydrides 
the  combustion  proceeds  through  an  initial  pyrolysis  reaction  which  is 
then  followed  by  the  combustion  of  the  metal  substrate.  The  result  is 
the  existence  of  multistage  reaction  sequence  in  which  the  first  stage 
has  a  relatively  cool  flame  temperature  which  is  then  followed  by  a  high 
temperature  secondary  combustion  stage.  Because  of  the  time  lag 
between  these  two  stages,  the  combustion  efficiency  of  metal  hydride 
particles  in  a  rocket  engine  is  lowered  relative  to  efficiencies  for  similar 
metal  particles  especially  for  short  engine  residence  times.  The 
extent  to  which  these  reactions  are  delayed  is  therefore  of  fundamental 
interest  to  propulsion  engineering. 

II.  EXPERIMENTAL  APPROACH 

Lithium  aluminum  hydride  -  LiAlH4  was  chosen  as  a  typical  light 
metal  hydride  for  the  study  of  ignition  and  combustion  processes.  This 
fuel  was  chosen  since  it  is  of  current  interest  for  some  propellant  sys¬ 
tems  and  secondly  because  its  thermal  degradation  and  physical  properties 
are  relatively  well  characterized.  Thermal  degradation  of  this  com¬ 
pound  proceeds  through  at  least  three  separate  steps.  1,2  Differential 
thermal  analysis  studies  shown  in  Figure  1  schematically  show  the 
proposed  reaction  sequence  for  this  degradation.  These  are  correlated 
with  experimental  differential  scanning  calorimeter  data  presented  in 
Figure  2.  These  experimental  results  clearly  show  the  stepwise  staging 
of  the  decomposition  reactions  which  occur  in  the  slow  pyrolysis  of  this 
metal  hydride. 

To  study  the  combustion  process  under  fast  reaction  conditions,  a 
series  of  experimental  spectroscopic  studies  were  performed  to  deter¬ 
mine  the  time  delay  between  the  ignition  of  the  first  stage  and  the  sub¬ 
sequent  reaction  sequence  for  fuel  oxidizer  mixtures  of  LiAlH4  -  Li2Oz 
and  LiAlH4-NH4N03.  Primary  experimental  variables  were  the  particle 
size  and  O/F  mixture  ratios.  The  experimental  procedure  consisted  of 
mixing  the  fuel  and  oxidizer  as  a  powder  in  the  appropriate  mixture 
ratio.  Ignition  of  the  "solid  propellant"  was  then  accomplished  by  either 
spark  or  laser  initiation.  The  reaction  sequence  was  followed  using  a 
fast  streak  spectrograph.  This  spectrographic  method  allows  a  detailed 
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Figure  2.  Differential  Scanning  Calorimetric  Analysis  for  LiAlH* 
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analysis  of  the  buildup  and  staging  of  the  reactions  important  to  the 
ignition,  pyrolysis  and  subsequent  combustion  of  these  powders. 

A  typical  series  of  spectrographic  records  for  the  combustion  of 
a  mixture  of  LiAlH4  -  Li202  is  given  in  Figure  3.  On  this  record  time 
"zero"  is  on  the  left  hand  side  and  reaction  proceeds  to  the  right.  The 
gross  features  of  this  record  include  a  single  Li  atomic  emission  line 
followed  by  a  wider  spectral  response  which  corresponds  to  the  first 
stage  reaction.  In  the  upper  region  a  careful  analysis  indicates  the 
probable  existence  of  self-absorption  bands  in  the  region  of  the  water 
bands.  The  second  stage  is  then  initiated  as  shown  by  the  large  increase 
in  the  blackbody  emission.  These  wider  bands  include  a  marked  black- 
body  continuum  together  with  a  reasonably  well-defined  self-absorption 
band  at  about  4800  A,  which  corresponds  to  strong  aluminum  oxide  bands. 
The  series  of  records  clearly  show  the  effect  of  particle  size  variation 
on  the  onset  of  these  processes.  The  particle  size  was  varied  to  provide 
a  semimonodisperse  system  having  a  particle  diameter  less  than  300, 

150  and  75  microns,  respectively.  As  the  particle  size  is  reduced,  the 
time  required  for  the  pyrolysis  reaction  is  reduced  from  time  zero  as  is 
the  time  delay  between  consecutive  stages. 

A  systematic  study  of  the  effect  of  particle  size  vs  time  delay  for 
these  reactions  was  performed.  Results  of  this  study  are  given  in 
Figures  4,  5,  and  6.  Figure  4  shows  the  variation  of  ignition  delay 
time  as  a  function  of  an  area  to  mass  ratio  normalized  to  a  75p  particle 
size.  Both  steps  of  the  reaction  are  shown  together  with  the  time  delay 
between  the  stages.  As  can  be  observed,  there  is  a  very  pronounced 
dependence  of  the  ignition  and  staging  delay  time  on  the  available  surface 
area.  Also  there  are  major  changes  in  the  shape  of  curvature  for  small 
particles  which  is  dependent  on  the  oxidizer  used.  This  latter  effect 
shows  the  strong  dependence  of  initiation  on  oxidative  environment 
although  insufficient  data  were  available  to  investigate  the  pyrolysis  of 
the  oxidizer.  Figures  5  and  6  show  the  effective  time  delay  to  ignition 
for  the  Stage  I  and  Stage  II  reactions,  respectively.  These  are  plotted 
as  a  function  of  O/F  mixture  ratios  for  the  LiAlH4  -  Li202  fuel  oxidizer 
combination.  Again  the  effect  of  particle  size  is  pronounced. 


Li202  (48)-LiAlH4  (48) 
O/F  =  3 


Li202  (48)-LiAlH4  (100) 
O/F  =  3 


Liz02  (100)-LiAlH4  (48) 
O/F  =  3 


Li202  (100)-LiAlH4  (100) 
O/F  =  3 


Li202  (200)-LiAlH4  (48) 
O/F  =  3 


Li2Oz  (200)-LiAlH4  (100) 
O/F  =  3 


TIME 


Figure  3.  Spcctrographic  Records  LiAlH4-Li202 
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Figure  3  (Cont) 


UAIH4  -  NH4NO3 
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IGNITION  DELAY  TIME  FOR  STAGE  I  REACTION 
LiAIH4  -  Li202 


Figure  5.  Time  Delay  to  Ignition  -  Stage  I  Reaction 


III.  THEORETICAL  APPROACH 


The  experimental  results  for  the  ignition  of  lithium  aluminum 
hydride  indicate  the  presence  of  two  stages  of  ignition.  These  two 
stages  plus  the  porosity  of  the  material  and  the  resulting  considerations 
of  coupled  mass  and  heat  transfer  within  the  material  must  be  included 
in  the  model  for  ignition.  Because  of  these  added  considerations, 
theoretical  models  describing  light  metal  hydride  ignition  must  be  much 
more  complex  than  the  current  models  for  the  ignition  of  molten 
metallic  particles. 


Based  on  the  experimental  results,  a  model  for  the  ignition  has 
been  formulate.*  and  is  presented  in  Figure  7.  The  main  steps  during 
the  transient  period  prior  to  final  ignition  include  the  following: 


Step  1:  Transient  Heating  -  The  particle  is  heated  up  from 

ambient  conditions  due  to  radiated  and  convective  heat 
transfer. 


Step  2:  Stage  I  Ignition  -  The  hydrogen  is  evolved  and  diffuses 
outward  from  the  particle  to  the  surface  where  ignition 
and  combustion  occurs. 


Step  3:  Lag  Time  -  After  the  completion  of  hydrogen  pyrolysis, 
a  is  provided  to  allow  time  for  the  oxygen  to  diffuse 
back  to  the  surface. 


Step  4:  Oxidation  -  As  the  oxygen  diffuses  to  the  surface,  the 
oxygen  reacts  with  the  exposed  aluminum  oxide. 

Step  5:  Vaporization  -  The  continued  heat  transfer  to  the 

interior  of  the  particle  from  the  environment,  and  the 
hydrogen  flame  plus  the  heat  generated  from  the  hydro¬ 
gen  evolution  and  aluminum  oxide  formation  result  in 
eventual  melting  and  vaporization  of  the  core. 

Step  6:  Stage  II  Ignition  -  The  vaporized  metal  diffuses  outward 
from  the  particle  and  ignites  when  the  proper  thermal 
balance  and  concentrations  of  aluminum  and  oxygen  are 
attained. 
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Flame 

Front 


Step  1.  Transient  Heating  Step  2.  Stage  I  Ignition 


Step  3.  Lag  Time 


Step  4.  Oxidation 


Step  5.  Vaporization  Step  6.  Stage  II  Ignition 

Figure  7.  Hypothetical  Model  for  Ignition  of  Light  Metal  Hydride 
Type  Particles 
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The  hypothetical  model  has  been  used  as  a  basis  for  a  digital 
computer  program  to  predict  the  Stage  I  and  Stage  II  ignition  times  and 
the  variation  of  the  hydrogen  flame  radius  with  time  during  hydrogen 
pyrolysis.  The  program  utilizes  the  nodal  approach  with  stepwise 
solution  of  differential  equations.  Reactions  external  to  and  within  the 
spherical  particle  are  assumed  to  be  symmetrical  so  that  the  nodes  are 
taken  as  concentric  "onion  skins"  within  the  particle.  Heat  and  mass 
balances  are  taken  across  each  node.  Whenever  reactions  occur  the 
changes  in  mass  and  energy  are  added  to  the  appropriate  node  within 
the  sphere.  The  hydrogen  pyrolysis  and  resulting  flame  zone  external 
to  the  particle  is  treated  by  an  appropriate  energy  balance  including  the 
enthalpies  of  the  reaction  and  the  reaching  and  product  gases.  The 
variation  of  the  radius  of  the  flame  front  is  predicted  by  assuming  a 
stable  flame  front  with  a  linear  variation  in  flame  radius  as  a  function 
of  the  total  enthalpy  difference  available  for  the  reaction.  According 
to  Brzustowski,  3  this  assumption  of  linear  variation  is  valid  under 
certain  conditions  for  stable  flame  zones. 


A  typical  heat  transfer  balance  across  a  node  within  the  particle 
is  illustrated  in  Figure  8.  In  this  figure,  the  conductive  heat  transfer, 
the  stored  energy  and  chemical  kinetic  energy  terms  are  all  included. 
The  energy  balance  becomes 


,  ^  ua  fdT\  UA 

(mCp>>  U72  1  137],-,  ‘ kAi 


(i) 


which  can  be  modified  by  making  the  assumption  that  the  derivatives 
can  be  approximated  by  the  ratio  of  small  increments  of  change 


This  assumption  allows  the  substitution  of  increments  into  (1)  to  obtain 


(t£  -  t2)  k 


[  Aj  (T i  -T2)  -  A2(T2  -T3)]  +  q. 


(4) 


which  can  readily  be  set  up  in  FORTRAN  IV  or  ALGOL  language  and 
solved. 

An  excellent  discussion  of  shell  balances  and  the  terms  for  heat 
and  mass  transfer  within  spheres  is  included  in  the  text  by  Bird, 

Stewart  and  Lightfoot.4  The  mass  diffusion  balance  using  node  concen¬ 
tration  levels  and  the  effective  diffusivity  form  is  completely  analogous 
with  the  heat  transfer  balance.  The  final  form  for  stepwise  differential 
solution  is  similar  to  (4) 

Cz  =  C2  +  (~j  [Al(Cl  -  C2)  -  A2(C2  -  C3)j 

(5) 

/  At  \  ,  n 

+  (-  J  kl 

where  kj  Cn  is  the  term  for  change  in  concentration  due  to  chemical 
reactions. 

For  the  heat  transfer  at  the  surface  of  the  particle,  a  quasi - 
steady -state  balance  is  made  between  the  heat  input  and  the  heat  con¬ 
ducted  into  the  sphere.  The  heat  input  includes  radiation  and  convection 
from  the  environment  and  radiation  from  the  flame  zone  when  hydrogen 
pyrolysis  is  occurring. 

For  the  mass  transfer  at  the  surface,  the  particle  is  assumed 
to  be  relatively  at  rest  with  respect  to  the  environment  so  that  forced 
mass  transfer  is  not  considered.  Although  the  mass  diffusion  to  and 
back  from  the  flame  'zone  is  recognized  as  a  major  factor  in  flame 
dynamics  and  ignition  delay  by  many  sources  including  Kuehl  and 
Zwillenberg,  5  the  external  diffusion  for  this  program  is  considered 
simply  as  the  diffusion  of  the  reacting  gas  outward  from  the  last  node 
to  the  surface.  Similarly,  for  diffusion  of  oxygen  inward  from  the 
surface,  the  oxygen  concentration  in  a  normal,  one  atmosphere  environ¬ 
ment  was  considered  as  the  driving  potential.  These  approximations 
should  not  significantly  affect  the  total  ignition  times.  However,  the 
time  to  ignition  for  Stage  I  may  be  predicted  slightly  prematurely  since 
no  time  delay  for  diffusion  of  the  hydrogen  into  the  surrounding  gaseous 
environment  is  included. 

The  computer  program  has  been  written  in  ALGOL  language 
since  an  extended  Burroughs  B5500  is  available  at  the  Denver  Research 
Institute.  This  program  has  been  named  as  the  PARTIG  Program  for 
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particle  ignition.  An  iterative  technique  was  used  for  solving  the  step¬ 
wise  differential  equations  (4)  and  (5)  for  heat  and  mass  transfer  at 
each  node.  Again,  the  heat  and  mass  transfer  reiterative  techniques 
are  analogous  so  that  the  computer  flow  diagrams  for  both  are  similar. 
The  reiterative  flow  diagram  for  heat  transfer  is  presented  in  Figure  9. 
The  iterative  procedure  involves  taking  heat  and  mass  balances  for  a 
particular  node  or  shell  and  then  progressing  to  the  next  node.  When 
all  nodes  are  completed,  the  present  time  for  the  values  of  temperature 
and  concentrations  for  all  nodes  is  switched  to  the  past  time  for  the 
next  time  step  and  the  nodal  iteration  procedure  is  again  repeated  for 
the  next  time  increment.  Note  that  although  the  explicit  solution  tech¬ 
nique  is  used,  oscillations  in  temperature  and  concentrations  are  pre¬ 
cluded  by  choice  of  a  sufficiently  small  time  increment  for  calculations. 
This  time  increment  is  estimated  by  the  time  constant  for  heat  transfer 
which  is 


where  all  terms  apply  to  a  particular  node. 


Program  Listing  and  Input  Data 


An  example  listing  for  the  input  data  and  equations  for  the 
PARTIG  Program  is  presented  in  Figure  10.  The  input  data  in  this 
listing  is  based  on  a  50p  lithium  aluminum  hydride  particle.  The  igni¬ 
tion  source  for  this  analysis  was  considered  as  the  hot  surrounding 
gases  in  a  solid  propellant  rocket  chamber  at  2500°K  since  the  potential 
use  of  these  particles  as  additives  in  solid  rocket  propellants  may  be 
desirable.  Also, for  comparisons  with  the  experimental  data  on  powder, 
the  reacting  temperatures  for  powder  and  the  hydrogen  flame  tempera¬ 
ture  would  be  close  to  2500 *K  (from  2200 'K  up)  and  would  contribute 
about  the  same  amount  toward  ignition  of  a  single  particle. 


For  the  hydrogen  evolution  and  heat  of  reaction,  the  experi¬ 
mental  data  of  Block  and  Gray1  as  presented  in  Figures  1  and  2  was 
used.  To  simplify  calculations,  the  hydrogen  evolved  and  the  heat  of 
reaction  were  considered  to  occur  as  steps  at  the  peak  temperatures 
and  the  total  amounts  of  hydrogen  and  heat  ev  lved  for  each  tempera¬ 
ture  region  indicated  in  these  figures  was  assumed  to  occur  at  the 
temperatures  of  the  peaks. 


Calculation  ot 
Flame  Enthalpy 
and  Radius 
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COMMENT  LI-AL-H4  PARTICLE  IGNITION  ; 

BEGIN 

FILE  OUT  LINE  4(2,15); 

A  M0063  DATELINE  ("EARTIG")  ; 

FORMAT  OUT  FMTl(E13.6.X4/ll(iO(F9.2,X2)/)); 

FORMAT  OUT  SEACEl( ) ; 

FORMAT  OUT  FMT2(E13.6,X4,I3,X4,E13.6,X4,E13.6,X4,E13.6)  ; 

FORMAT  OUT  FMT3(E13-6,X4,E13.6,X4,E13.6); 

FORMAT  OUT  FMT4(X10, "SHELL  TEMPERATURES"); 

FORMAT  OUT  FMT5(X10, "OXYGEN  CONCENTRATION" ) ; 

FORMAT  OUT  FMTd(X  10, "ALUMINUM  VAPOR  CONCENTRATION"); 

FORMAT  OUT  FMT7(X10, "HYDROGEN  PYROLYSIS"); 

FORMAT  OUT  FKr8(X10,"LIALH4  PARTICLE  IGNITION,  INITIAL  VALUES"); 

FORMAT  OUT  FKT9(E13.6,X4/4(6(E13.6,X2)/)); 

FORMAT  OUT  FMT10(X10, "ALUMINUM  OXIDE  EQUIVALENT  <f2  COKCEN."); 

INTEGER  J,I,M,I1,J1,B; 

INTEGER  Y,Z,D; 

REAL  DELE,DELT,R,K1A,RH01,CP1A,K2A,RH02,CP2A,K3A,RK03-CP3A, 

PI,TIME,TIME2 ,COHD,ALFA,CP,RHO ,TF1,TF2 ,TF3 ,TPGI ,TPG2 , 
TPG3,H.DELH1,G,H2,!»,H01,CI,CW,RF,HP,FEPG,FES,FE1,FE2, 

FE3  ,SIGMA,HPG,HjJ25  ,TPG,QIAS  ,CPH; 

REAL  P0,LAG,CAC1,CAC2,CBC1,DIFA,DIFB,MELT,LV,CAD,CA0,AI/,MAL; 

REAL  MO; 

ARRAY  RADC0:21,0:2lJ,  AC0:2  ,0:^J  ,DV^):2i]  ,T[p:2,0:22j  ,HTlp:2]  ; 

ARRAY  DELH(&:2,0:22],HDB>:2,0:2g  ,AL[(0:2,0:22|  ,ALTQ):2,0:23  ; 

ARRAY  CAO):2,0:223  ,C3 [0:3,0:222  ,TIMQ):2]  ,Q0):^  ,vl<0:S| ; 

ARRAY  /C0:2, 0:223; 

LABEL  PI; 

IABEL  P2  ; 

FOR  J«-l  STEP  1  UNTIL  6  DO 

BEGIN  OCljjl^-O; 

CBCl,fl-*-0.0; 

caD>jj+-o; 

TCL,^J-*-300; 

TIME«-0.0 

END; 

WRITE(LINE,FMT8) ; 

WRITE ( LINE, FMT4) ; 

WRITE(LINE,FMT1, TIKE, FfR  J-*-l  STEP  1  UNTIL  6  DO  tQ.,j]  ); 


Figure  10.  Listing  for  the  PARTIG  Program  -  As  Run  for 
LiAlH4  Ignition 
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WRITE 'LINE, FMT5); 

WRITE(LINE,FMT9,TIME,F2)R  J-l  STEP  1  UNTIL  6  M 

caD-.J]  ); 

WRITE ( LINE , FMT6 ) ; 

WRITE(LINE,FMT9,TIME,F@R  J— 1  STEP  1  UNTIL  6  Ut 

cbCi.J]); 

WRITE(lINE,SPACE1); 

DELE  4*0. 200; 

DELT«-1.00@-6  ; 

R*-2.5@-3; 

Y4-0; 

Z4-0; 

DIF.44-  0.200; 

DIFB4-  0.150; 

Vt  4-1.0; 

LAG4-25 ; 

CAC14-0; 

CAC24-7.8l@-10; 

CBC14-1.04@-S; 

MELT4-991; 

M0 4-16.0; 

Q[I3^0; 

CA04-8.2@-7; 

Vp40; 

AI48 4-3900; 

LV 4-2360 ; 

MAL4-26.97; 

D4-j.; 

KIA4-9.09§-3; 

RH/14-0.917; 

CP1A4  0.200; 

K2A4-9.09®-3; 

RH0 24- 0.820; 

CP2A4-  9 . 200 ; 

SIGMA4- 1.365®- 12; 

PI4-3.1U  16; 

TPG14-2500; 

TPG24-0; 

TFG34-0; 

TF14-  443 ; 

TF24-2320; 

TF34-0  ; 

CAD4-RH02»MO/(2  .C*!AL) ; 

H0254-O  ; 

CP04. 0.262  ; 

FEPG4-1.00  ; 

FES4-  C .  90  5 
CPH4-3.25  ; 

HP4-1.43S4; 

14-2; 

114- I-l; 

TFG4-TRJ1  ; 

HPG4*1.6@-2  ; 

QLAS  «-0.0; 

F«!R  J4-1  STEP  1  UNTIL  6  Vf 


Figure  10.  (Continued) 
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BEGIN 

Jl+J-l; 

RAD[J1,1]*-R8(l-(J1-1)*DELE); 

RAD  &T1 ,2  R®(  1-J10DELE  ) ; 

A[jl,l]*-li«PlOR.\DtJl,q*2  ; 

ACJl,2j*-4«Pl«RADtJl,23*2  ; 

DVCJll*-  4/3*(RADCJl,ll*3-RADCri,22*3)#n; 

IF  J=2  THEN  B*-2  ELSE  B*-l  ; 

IF  TCr.l,J3  =  723  THEN 

BEGIN 

C«)ND*-K2A; 

CP-*-CP2A; 

RH0*-RH02; 

END  ELSE 
BEGIN 

C0ND-*-KlA; 

CP-4-CP1A; 

RH0*-RH01; 

END; 

if  tDci,j^tfi  THEN 
BEGIN  DELH[l,J3^*56.9>  Cl«.1.88®  5 
END; 

IF  TCll,J3i483  THEN 
BEGIN  DELHfl.JD '^■115.4;  CL *-8.51®  4 
END; 

if  t&i,«0*548  then 
BEGIN  DELH[I,J3^166.6;  Cl**'4.86®  4 
END; 

IF  T[ll,J]^723  THEN 
BEGIN  DELHCr,J]-*-327.6;  Cl-*-1.45@  4 
END; 

IF  J=1  THEN  H**0  ELSE 

BEGIN 

H‘*-DVCJ3®(DELHflJjJ  -DELH[ll,J3  )  5 

H2<*-CI»H  ; 

Hj5*-CP0®(TK)-298)  +H025; 

HD[I,J]-*- H2»{  HP-0 . 5K»-CPMj(  T  ^1  ,Jl -298) ) ; 

IF  H*0  THEN 

BEGIN  WRITE { LINE, FMT7); 

WRITE ( LINE , FMT2 ,TIME,J1,H2 ,HDC[ ,  J} ) ; 

END; 

END; 

IF  T&1,61>723  THEN 

BEGIN 

TIK  DAYTIME; 

Y^-l; 

IF  TIME  £  TIM  [03+IAGQDELT  THEN 
BEGIN  Z«-1.0; 

if  caCci,j3^caci  THEN 
BEGIN  Q[I3*-1.0; 

IH-2; 

<>*-2®(QCl]  -Q&lJ  )»DV fcl)^L0*M0®( CAD- j) £,,£)> 

END;  IF  J=1  THEN  CA[i,J3*-2®CA0 

EISE  CA[I,J]-** Z»{CAD;i,J]  +DIFA»DELT/( DELEtDVCJll )  (ACll.13 
®(CACll,Jl]-CAtLl,J])-ADl,2]«(CADl,0Q-CAGl,J+i3))/R) 


Figure  10.  (Continued) 
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-(Q  &]-Q  Cri]  )®(CAD-^  C2,  J3) ; 

IF  CA[I,J3  <0  THEN 

BEGIN  +Zfc(CA til ,  J]  +DIFA«DELT/( DELEttDV tjl]« R ) 

«(AUJ1 , 1]©(CA  fll  -CA  CI1 ,  JU ) 

-ADl ,2 38(CALH , J] -CAQ1 , J+13) ) ) ; 

IF  <J[1,J3>CAD  THEN  CA[l,J>  0C1, JO-CAD  EISE 
BEGIN  CAri.J^^O; 

QtI3^*0; 

G  «*-2*DV  CJll«AW».f^8(0  C2,j]  -0D-,J3) ; 

END;  0[l,J]-*-0C2,J3; 

END; 

END; 

END; 

IF  T(Tl,J3  5-MELT  AND  D=1  THEN 
BEGIN  VCl]-*-1.0; 

IF  ALTTIIjJ}*  LV»DVCJll©RH02  THEN  AL[T,j3-*-0  ELSE 
BEGIN  IF  J=1  THEN  AL[l,Jj-*-0  ELSE 

ALfl  C0ND/(  RSOELE  )€)(B8A  [J1 , 118(  Till  ,Jll -T  Til  ,J]  ) 

-A  [J1 ,2l8(  T[ll,Jl-TCll  ,  J+Il ) )®  DELT ; 

END; 

CB[I  ,J] ■*-  CB  [II  .J]  - DIFA9DELT/ ( DELEtDV  Dll  )3(A GJ1, ll 
®(CBCIl,Jn-CBtll,J])-ACJl,236(CBCll,J]-CBCIl,J+13))/R 
+(Vtn-VCIll)®l'iAlAALtI,J3/(DVCJl'y*LV) ; 

ALTtI ,  J3*-ALTLI1 ,  J]  +ALCI ,  JJ ; 

END  ELSE  AL[I,J1*0; 

ALFA"*- C0ND/( RH00CP)  ; 

IF  J=1  THEN 

T[I,J3**-  TCll,2l  +DELE/(  28C®ND  )&(HPS®(  TPG-T  til .  13  ) 

+SIGM/'®(  FEB3©rPG*l*-FES0TtIl ,  13*4+FE1£TF2*4 )  )8R 
ELSE 

TtI,J>-  T[I1,J]  +ALFA3DELT / ( R0DELE&DV tj  13  )  ®  ( BtA [ J 1 , lT® 
(Ttll,Jl3  -TCll , J3)-ACJ1,23®(TCE1, J3 -TCll, J+tf ) ) 

+DEL?/(  RH4®CPS0V[J13  )®(H+G-ALCl,  J]) ; 

END; 

TIME2-*-TIME  +  DELT  ; 

HTm-HDtl  ,11+HDCl  ,21 +HDCI  ,31 +HDCI  ,4]+HDCI,53+HDlI  ,63; 

IF  HTtn-HTtIll=0  THEN  RF«- P.  ELSE 

RF*  RaHTtll  /( HT  til -HTCI13  jSiyf ; 

WRITE  ( LINE , FMT3 ,TIME2 ,HT III , RF ) ; 

WRITE ( LINE , FMT4 ) ; 

WRITE ( LINE , FMT1.TIME2 , FGR  J-v-1  STEP  1  UNTIL  6  D 0  TtI,J3); 
WRITE (LINE, FMT5) ; 

WRITE ( LINE, FMT9,TIME2  ,F0R  J«-l  STEP  1  UNTIL  6  D0 
CALI,J3); 

WRITE ( LINE , FMT 10 ) ; 

WRITS(LINE,FM?9,TIME2,F0R  J*-l  STEP  1  UNTIL  6  D0 
0tl,J3); 

WRITE ( LINE, FMT6) ; 

WRITE ( LINE , FMT9 ,TIME2 ,  PGR  J  1  STEP  1  UNTIL  6  DG 
CBCI,J3); 

WRITE ( LINE , SPACE1 ) 

PGR  M*-l  STEP  1  UNTIL  6  DG 
BEGIN  Tll,M3«-TL2,Ml ; 

HT[13*-HTC21 ; 


Figure  10.  (Continued) 
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P2: 


DELH  [1  DELHte  ,Kl ; 

TIME*-TIME2; 

CA  Cl  ,M3«-CA£2  jMl ; 

CBtl,M]*-CBC2,M3; 

ALTU,Ml*-ALTt2,M}; 

1IDC1,M]4-HDC2,M1; 

Qtl3^-QC23 ; 

V[11^VC23; 

END; 

IB’  TIME  DELT^  10000  THEN  G0  It  PI; 
DATELINE (0)  ; 

END. 


Figure  10.  (Continued) 
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For  the  program,  the  input  values  selected  for  thermal  and 
physical  properties  are  indicated  on  the  listing  (Figure  10)  on  cards  36 
through  78.  The  effective  mass  diffusivities,  DIFA  and  DIFB,  were 
estimated  from  values  for  molecular  diffusivity  tabulated  in  Bird, 
Stewart  and  Lightfoot.  4  These  values  for  molecular  diffusivity  cor¬ 
rected  to  effective  diffusivity  for  the  porous  pellet  by  using  effectiveness 

6  3/2 

corrections  from  Petersen.  Also,  a  temperature  correction  of  T  f 

for  the  average  temperature  was  corrected  to  effective  diffusivity  for 

the  porous  pellet  by  using  effectiveness  corrections  from  Petersen.  6 

Also,  a  temperature  correction  of  T3/2  for  the  average  temperature 

was  applied  based  on  the  Chapman-Enskog  equation  as  presented  in 

Bird,  Stewart  and  Lightfoot4  assuming  a  perfect  gas. 

The  critical  concentrations  were  not  available  for  a  2500 °K 
environment  and  some  of  the  other  properties  were  estimated.  Thermal 
and  physical  constants  such  as  heats  of  reaction  and  vaporization, 
thermal  conductivity,  and  density  were  based  on  data  from  Perry,  7  the 
CRC  Handbook8  and  the  JANAF  Tables.  9  Since,  during  the  pyrolysis 
of  hydrogen  and  other  steps,  it  is  not  clear  whether  or  not  the  total 
lithium  content  is  pyrolyzed,  the  melting  temperature  for  the  core, 
MELT,  is  based  on  two-phase  data  on  the  melting  temperature  of  Li-Al 
from  Hansen.  10 

The  convective  heat  transfer  coefficient  HPG  was  estimated  for 
free  convection  assuming  that  the  particle  was  at  rest  relative  to  the 
surrounding  environment.  According  to  McAdams,  11  data  has  been 
presented  which  indicates  that  for  small  particles  (below  ~  l/4  inch  in 
diameter)  the  Nusselt  number  approaches  a  constant  value  which  is 
close  to  or  actually  the  transfer  value  due  to  gaseous  conduction. 
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IV.  DISC  JSSION 

The  results  of  this  investigation  are  in  general  agreement  with 
the  kinetic  mechanism  proposed  by  previous  investigators.1’2  Garner 
and  Haycock  investigated  the  isothermal  decomposition  of  Lithium 
Tetrahydraluminate  and  found  a  reaction  sequence  corresponding  to  the 
formal  equation 

1 .  LiAlH4  -  Li A1H2  +  H2 

2.  LiAlH2  -  LiH  +  A1  +  |H2 

Block  and  Gray  modified  this  reaction  sequence  by  including 

3.  LiH  -  Li  +  £H2. 

These  workers  utilized  a  slow  pyrolysis  during  which  the  sample  was 
heated  at  a  rate  of  temperature  rise  approximately  10°  1  min.  Although 
the  present  work  has  a  much  higher  thermal  rise  rate,  the  general 
reaction  sequence  appears  to  be  basically  similar.  The  spectrographic 
studies  are  characterized  by  the  onset  of  an  atomic  lithium  emission 
line  at  relatively  low  temperatures.  This  emission  remains  at  a  con¬ 
stant  "temperature"  for  a  considerable  length  of  time.  There  then 
follows  a  discrete  increase  in  the  intensity  of  emission  in  which  addi¬ 
tional  atomic  lithium  lines  are  activated  indicating  an  increased  reaction 
temperature.  This  condition  forms  a  discrete  second  stage  which  is 
then  followed  by  an  afterburning  of  the  residual  substrate.  The  sequence 
supports  a  reaction  mechanism  which  can  be  described  by  the  following 
equations: 


la. 

LiAlH4  - 

LiAlH2  +  H2 

Step  1  <  lb. 

LiAlH4  - 

LiH  +  AIH3 

__  lc. 

LiH 

Li  +  |H2 

Step  2  /”  ^a‘ 

LiAlH2  - 

LiH  +  |H,  +  Al 

2b. 

LiH 

Li  +  H2 

Step  3 

A1  +  Oz  - 

ai2o3  . 

The  principle  difference  in  this  mechanism  is  that  it  is  necessary 
to  include  an  additional  early  step  which  can  produce  atomic  lithium. 
This  could  entail  a  diffusion  of  the  lithium  atom  from  an  internal  site 
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to  the  surface  or  a  disproportionation  to  form  lithium  hydride  and  alu¬ 
minum  hydride.  Alternatively  the  early  intense  lithium  emission  could 
also  be  explained  as  "steady  state"  equilibrization  of  reactions  1  and  2 
of  Block  and  Gray.  However,  if  this  were  occurring,  then  a  more 
gradual  buildup  to  the  second  stage  of  the  reaction  would  be  anticipated. 

Results  for  the  PARTIG  Program  are  shown  in  Figure  1 1  in 
which  time  and  temperatures  and  concentrations  are  indicated  for  each 
node.  Also,  the  hydrogen  evolved  and  heat  of  evolution  are  listed  at  the 
appropriate  peak  temperature  in  Figure  12.  Due  to  the  use  of  discrete 
peak  temperatures  the  hydrogen  evolution  appears  as  "spikes."  An 
estimate  of  the  flame  shape  for  a  more  continuous  modified  Arrhenius 
hydrogen  evolution  is  indicated  by  the  dashed  line  on  the  figure.  The 
Phase  I  ignition  would  appear  to  start  at  about  0.  3  milliseconds. 

Theoretical  prediction  for  the  temperatures  of  the  surface  and 
at  the  center  of  the  particle  are  presented  in  Figure  12.  The  results 
indicate  a  linear  rise  in  temperature  with  time  with  no  visible  perturba¬ 
tions  (at  the  scale  of  the  figure)  of  temperature  due  to  heat  generated 
from  chemical  reactions  within  the  nodes.  It  is  noted  that  the  lag 
between  the  temperature  of  the  surface  and  the  center  of  the  particle 
is  almost  constant  at  only  10oK.  The  lag  time  and  diffusion  of  oxygen 
into  the  surface  after  hydrogen  pyrolysis  did  not  appear  to  affect  the 
temperature.  Also,  since  the  concentration  of  oxygen  in  the  atmosphere 
is  rather  low,  the  formation  of  aluminum  oxide  was  limited  to  an 
extremely  thin  layer.  Because  the  oxygen  diffused  into  the  particle  was 
consumed  rapidly  near  the  surface,  no  concentration  was  available  for 
reaction  with  the  metal  vapor  as  the  core  reached  melting  temperature. 
Therefore,  the  Stage  II  ignition  was  assumed  to  take  place  as  soon  as 
the  vaporized  metal  began  to  diffuse  to  the  surface.  Although  the  991  °K 
melting  temperature  is  probably  not  sufficient  for  ignition,  the  2500 “K 
surrounding  environment  was  considered  adequate  to  almost  instan¬ 
taneously  bring  the  metal  vapor  to  the  ignition  temperature.  Based  on 
this  assumption,  the  Stage  II  ignition  occurs  at  about  1.6  milliseconds. 

The  reduced  experimental  results  as  presented  in  Figures  4, 

5  and  6  were  compared  with  the  analytical  predictions  (Figure  12  and 
13)  and  the  results  are  presented  on  Figure  14.  The  50p  particle  delay 
times  from  the  program  were  normalized  to  the  75p  diameter  by  use  of 
the  D2  law  for  comparison  purposes.  The  comparison  indicates  that  the 
analytically  predicted  Stage  I  and  Stage  II  ignition  times  are  in  close 
agreement  with  the  observed  experimental  data  for  an  O/F  ratio  of  1.0. 
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IIAIH*  PARTICLE  IGNITION#  iNITtAL  VALUES 
SHELL  TEMPERATURES 
O.OOOOOOMOO 

300,00  300,00  300.00  300.00  300.00  300,00 

OXYGEN  CONCENTRATION 
0.0000009*00 

0.0000009*00  0,0000009*00  0.000000**00  0.0000009*00  0.0000009*00  0,0000009*00 

aluminum  vapor  concentration 

0.0000009*00 

0.0000009*00  0,0000009*00  0,000000**00  0.000000**00  0 • 000000**00  0,0000009*00 

4,0000009*07  0.0000009*00  2.5000009*03 

SHELL  TEMPERATURES 
4.0000009*07 

302.43  300,00  non, oo  300.00  300.00  300,00 

OXYGEN  CONCENTRATION 
4.0000009*07 

0.000000**00  0,000000**00  0.000000**00  0.000000**00  0.0000009*00  0.0000009*00 

ALUMINUM  vapor  CONCENTRATION 
4.0000009*07 

0.000000**00  0,000000**00  0.0000009*00  0.000000**00  0.000000**00  0.0000009*00 

a, 0000009*07  0.0000009*00  2,5000009*03 

SHELL  TEMPERATURES 
8.0000009*07 

302.43  301.49  300,00  300.00  300.00  300.00 

OXYGEN  CONCENTRATION 
8.0000009*07 

0.0000009*00  0.0000009*00  0.0000009*00  0.000000**00  0.0000009*00  0,0000009*00 

aluminum  vapor  concentration 

8.0000009-07 

0.0000009*00  0.0000009*00  0.0000009*00  0.000000**00  0.0000009*00  0,0000009*00 

l,2000009*0«  0.0000009*00  2,5000009*03 

SHELL  TEMPERATURES 
1.2000009-06 

303.92  301,78  300.48  300.00  300,00  300.00 

OXYGEN  CONCENTRATION 
1,2000009*06 

0,000000**00  0,000000**00  0.0000009*00  0.000000**00  0,0000009*00  0,0000009*00 

aluminum  VAPOR  CONCENTRATION 
1.2000009-06 

0.0000009*00  0,000000**00  0,000000**00  0.000000**00  0,0000009*00  0.0000009*00 


1 ,6000009*06  0.0000009*00  2.5000009*03 

SHELL  TEMPERATURES 
• ,6000009*06 


Figure  11.  PARTIG  Program  Printout  Results  -  Example  for  LiAlH4 


H2  I.ibi 


Mtli] 


-  Particle  Diameter:  50{i 

-  Environment  =  2500*K,  1  ATM 

-  No  Forced  Convection 
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Although  this  comparison  appears  to  confirm  the  delays  estimated 
by  the  PARTIG  Program,  the  analytical  and  experimental  conditions 
differed  considerably  for  oxidizer.  Because  of  these  differences,  addi¬ 
tional  analytical  solutions  with  the  PARTIG  Program  will  be  needed  in 
order  to  completely  verify  the  accuracy  of  the  program  for  predicting 
ignition  times. 
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V.  CONCLUSIONS 

A  number  of  conclusions  can  be  formulated  even  with  limited 
data  presently  available. 

1.  This  simple  combustion  model  used  to  describe  the  com¬ 
bustion  of  molten  metal  spheres  requires  modification  to 
account  for  a  multistage  reaction  sequence  in  which  the 
"cool  flame"  pyrolysis  is  followed  by  a  "hot  flame"  com¬ 
bustion  zone. 

2.  The  time  lag  required  for  the  ignition  of  the  particle  is 
strongly  dependent  on  the  particle  size  and  surface  aree.s 
and  varies  as  a  function  of  the  O/F  mixture  ratio. 

3.  Prediction  of  ignition  delay  times  for  multistaged  ignition 
in  porous  light  metal  hydride  particles  can  be  described 
using  a  nodal  shell  balance  analysis. 

4.  The  analytical  and  experimental  results  are  in  reasonable 
agreement  and  tend  to  confirm  the  assumed  theoretical 
model  for  multistage  pyrolysis  and  combustion  of  light 
metal  hydride  particles. 

From  an  engineering  point  of  view,  the  results  of  this  program 
indicate  that  the  stepwise  pyrolysis  of  a  metal  hydride  must  be  considered 
in  order  to  insure  efficient  combustion.  For  example,  the  optimization 
of  a  propellant  surface  regression  rate  is  dependent,  in  part,  on  the 
time  required  for  the  pyrolysis  of  the  hydride  component.  The  reported 
information  can  be  used  in  two  ways:  (1)  to  determine  the  maximum 
particle  size  necessary  for  the  complete  combustion  of  the  slow-burning 
fuel  component,  and  (2)  to  control  ignition  spiking.  The  necessary 
particle  size  can  be  estimated  for  these  propellants  from  an  examina¬ 
tion  of  Figure  4.  Ignition  spiking  can  be  reduced  by  placing  a  hydride 
layer  on  the  propellant  surface  and  using  the  built-in  reaction  stages  to 
control  pressure  rise  rates  at  a  minimal  impulse  loss.  An  estimate  of 
the  time  delay  can  be  developed  through  reference  to  Figures  5  and  6. 
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APPENDIX  A.  NOMENCLATURE 


A.  Text 


A 

C 

C' 

CP 

De 

k 

kj 

m 

n 

o/f 

q 

R 

Ri,  R2  •  •  • 
t 

tc 

T 

t' 

AV 

x 

P 

P 


Area 

Molar  concentration 

Molar  concentration,  present  time  interval 

Specific  heat  at  constant  pressure 

Effective  mass  diffusivity 

Thermal  conductivity 

Constant  for  reaction  kinetics 

Mass 

Order  of  reaction 
Oxygen/fuel  ratio 
Heating  rate  from  reactions 
Radius  of  particle 
Radii  of  nodal  surfaces 
Time 

Time  constant  for  conductive  heat  transfer 
Temperature 

Temperature,  present  time  interval 
Incremental  volume  of  node 
Distance  between  node  centers 
Density 

Ratio  of  nodal  shell  rings  to  total  radius,  AR/R 
Micron 


Subscripts: 


1,2.  .  .  Nodes 
melt.  Melting 


B. _ PARTIG  Program  Nomenclature 


A[J1, 1] 
A[J1, 2] 
ALFA 
AL[I,  J] 
ALO 

ALT[I,  J] 


Surface  area  into  shell  J  -  1 
Surface  area  from  shell  J  -  1 
Thermal  diffusivity 

Incremental  amount  of  metal  vaporized 
Heat  of  formation  of  aluminum  oxide 
Total  amount  of  metal  vaporized 


B 


Integer,  B  =  2  for  J  =  1,  E  =  1  otherwise 


Cl 

ca[i,  j  ] 

CAC1 
CAC2 
CAD 
CB[l,  J] 
CBC 1 
COND 
CPI  A 
CPI  B  S 
CP1C  J 
CP2A  1 
CP2B  > 

CP2C  J 

CP3A 
CP3B  > 
CP3C  J 
CPH 
CPO 

D 

DELE 

DELHI 

DELH[l, 

DELT 

DIFA 

DIFB 

DV[J1] 

FE1 

FE2 

FEPG 

FES 

G 

H 
H2 
HO 
H025 
HD[l,  J] 


Ratio  (H2/A1H,) 

Molecular  concentration  of  species  A 
Critical  concentration  of  A  for  first  ignition 
Critical  concentration  of  A  for  second  ignition 
Concentration  capacity  of  reaction  with  A 
Molecular  concentration  of  species  B 
Critical  concentration  of  B  for  first  ignition 
Thermal  conductivity 

Specific  heat  constants  for  material  #1  with 
CP  =  CPI  A  +  CP1B(T)  r  CP1C(T)2 

Specific  heat  constants  for  material  *2  with 
CP  =  CP2A  +  CP2B(T)  +  CP2C(T)2 

Specific  heat  constants  for  material  #3  with 
CP  -  CP3A  +  CP3B(T)  +  CP3C(T}3 

Specific  heat  of  hydrogen 
Specific  heat  of  oxygen 

Integer,  D  ^  2  for  Al2Oj,  D  =  1  otherwise 
Incremental  distance  between  shells 
Enthalpy  change  for  evolution  of  H2,  cal/gm 
Enthalpy  change  for  Jth  node,  cal/gm 
Incremental  step  change  in  time 
Mass  diffusivity  of  species  A 
Mass  diffusivity  of  species  B 
Incremental  volume  for  J1  shell 

Combined  emissivity  factor  for  H2  flame  front 
Combined  emissivity  factor  for  A1  flame  front 
Combined  emissivity  factor  for  propellant  gases 
Combined  emissivity  factor  for  particle  surface 

Enthalpy  change  for  oxide  formation 

Enthalpy  change  for  H2  pyrolysis 
Amount  of  hydrogen  liberated 
Enthalpy  of  02  for  H2  pyrolysis 
Base  value  of  Oz  enthalpy  at  25 °C 
Enthalpy  of  products  for  II2-02  reaction 


HP 

HPG 


H2  pyrolysis  reaction  energy 

Heat  transfer  coefficient  for  propellant  gases  to  particle 
surface 

HT[l]  Total  enthalpy  for  H2-O2  reaction 


I  Integer  for  present  time  increment 

II  Integer  (I  -  1)  for  previous  time  increment 

J  Integer  for  j'**1  particle  shell 

J1  Integer  for  J  -  1  particle  shell 


K1A 
K1B  - 
K1C  j 
K2A  'j 
K2B  S 
K2C  J 

K3A  1 
K3B  \ 
K3C  j 


Thermal  conductivity  constants  for  material  #1 
COND  =  K1A  +  K1B(T)  +  K1C(T)2 


Thermal  conductivity  constants  for  material  #2 


Thermal  conductivity  constants  for  material  #3 


LAG  Number  of  lag  time  increments  before  oxygen  back -diffusion 

LV  Latent  heat  of  vaporization  of  metal 

M  Integer  for  time  conversion  from  I  to  II 

MAL  Molecular  weight  of  aluminum 

MELT  Melting  temperature  of  metal 

MO  Molecular  weight  of  oxygen 

PI,  P2  Labels 

PI  pi  (3.  1416) 

PO  Partial  pressure  ratio  of  oxygen  to  1  atmosphere 

Q[l]  Switching  integer  {1  or  0) 

QLAS  Laser  heat  input 

R  Radius  of  particle 

RAD[J1,  1]  Radius  into  J1  shell 

RAD[J1,2]  Radius  inwards  from  J1  shell 

RF  Flame  front  radius  for  H2-O2  reaction 

RHO  Density 

RHOl  Density  of  material  #1 

RH02  Density  of  material  #2 

RH03  Density  of  material  #3 


SIGMA 


Stefan-Boltzmann  constant 


T[I,  J] 
T[I1,  J] 
TIM[I] 
TIME 
TIME  2 
TF1 
TF2 
TF3 
TPG 
TPG1  j 
TPG2  J 
TPG3  J 

V[I] 

Y 

Z 


Note: 


Temperature  of  the  shell  at  the  present  time 
Temperature  of  the  jth  shell  at  the  previous  time 
Time  at  end  of  hydrogen  pyrolysis 
Present  time 

Next  time  (TIME  +  DELT) 

Temperature  of  reaction  for  H2  pyrolysis 
Temperature  of  reaction  for  oxide  formation 
Reaction  temperature  for  metal  ignition 
Environment  gas  temperature 

Constants  for  propellant  gases  temperature 
TPG  =  TPG1  +  TPG2  (TIME)  +  TPG3  (TIME)2 


Switching  integer  (1  or  0) 
Integer 

Switching  integer  (1  or  0) 


Any  consistent  system  of  units  may  be  used  for  the 
program,  c  .  g.  s.  units  were  used  by  the  authors. 
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— T-hc  slow  pyrolysis  of  lithium  aluminum  hydride  proceeds  by  a  successive 
dehydrogenation  of  the  fuel  in  three  stages  over  a  temperature  range  between* 

.  210*  and  460®G.  The  combustion. of . lithium, aluminum  hydride  follows  a.  Simula  t* 
multistage  process  in, which  the,  pyrolytic  dehydrogenation  is  followed  by  an  after¬ 
burning  of  the  metal  substrate.  Experimental  studies  were  performed  -  to-  deter? 
mine  the  delay  times  between  successive  stages,  for  solid  mixtures  of  lithium 
.aluminum  hydride  with  selected  oxidizers.  Time  resolved  spectrographic 
-techniques  were  used  to  determine  the  duration  of  these  stages.  Studies  v/ere 
performed  using  Li2d2  and  NH4X03  oxidizers  over  a  particle  size  range  between 
75-300  microns  and  O/F  mixture  ratios  of  between  0.  5  and  3.  0.  Ignit'on  delay 
times  for  the  initial  reaction  were  found  to  vary  between  5  and  40  milliseconds  for 
the  first  stage  and  10  to  120  milliseconds  for  the  second  stages.  An  analytical 
model  was  developed  to  describe  the  observed  multistage  combustion  prccess. 

The  model  anc  .resulting  PARTIG  computer  program  utilize  the  nodal  approach 
with  stepwise  solution  of  the  equations  for  heat  and  mass  transfer/  ''Ijie._ 
analytically  predicted  ignition  delay  times  for  a  50p  particle  agreed  ilavcvgari.y 
with  the  experimentally  obtained  times.  *  , 
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